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Abstract 

We show that new universahty of Lyapunov spectra {Aj} exists in 
Hamiltonian systems with many degrees of freedom. The universality appears 
in systems which are neither nearly integrable nor fully chaotic, and it is 
different from the one which is obtained in fully chaotic systems on 1- 
dimensional chains as follows. One is that the universality is found in a 
finite range of large i/N rather than the whole range, where N is degrees of 
freedom. Another is that Lyapunov spectra are not straight, while fully chaotic 
systems give straight Lyapunov spectra even on the 3-dimensional simple cubic 
lattice. The universality appears when quadratic terms of a potential function 
dominate higher terms, harmonic motions are hence regarded as the base of 
global motions. 

1. Introduction 

In the last a few decades Hamiltonian systems with many degrees of freedom 
have been numerically investigated by integrating equations of motions. 
Hamiltonian systems being the foundation of statistical mechanics, one 
direction of the investigations is to check ergodic property originating in Fermi- 
Pasta- Ulam problem Another direction is to study their dynamical 
properties, which are transition from nearly integrable systems to stronger 
chaotic ones dynamical properties of phase transition structure of 
phase spaces [§] and so on. 

Here we focus on dynamical properties, in particular, universal structures 
of phase spaces which are not affected by details of systems. One of the 
structures is self-similar structure which is based on Poincare-Birkhoff theorem 
[0 in nearly integrable systems. Although this theorem is available in systems 
with two degrees of freedom, the self-similar structure is also supposed in 
systems with many degrees of freedom. 
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Lyapunov spectrum is usually used to study instability along a sample 
orbit. Moreover, it is useful to study structure of phase spaces both in 
dissipative |10[|Tl| and Hamiltonian ||5[||^ systems with many degrees of 
freedom, since it includes information on all directions in phase space. 
Although Lyapunov spectra in a system reveal dynamical properties of the 
system, we are interested further in properties which are not affected by details 
of systems. To detect such properties a useful approach is to find the universal 
form of Lyapunov spectra which is obtained in all the systems and whose cause 
indicates the properties. A universal form of Lyapunov spectra L{i/N) = Aj 
is reported in Hamiltonian systems which are 1-dimensional chains consisting 
of nonlinear oscillators [|12], where N is degrees of freedom. This universality 
gives the straight form for Lyapunov spectra, namely L{i/N) = a + b-i/N, and 
the straight form is also obtained with random matrices |jl2|] fl^ . Consequently, 
systems having straight Lyapunov spectra are regarded as fully chaotic ones 
which have not structure in phase spaces. 

In this paper, we show new universality of Lyapunov spectra appearing in 
systems with moderate strength of chaos which are neither fully chaotic nor 
nearly integrable, and study the cause of the universality. Here we define the 
word "universality of Lyapunov spectra" as L{i/N) approximately takes the 
same form in a finite range of large i/N regardless of total energy and details 
of systems. Systems with moderate strength of chaos are interesting for the 
following two reasons: one is that they have some unsolved problems and show 
interesting phenomena, for instance, second order phase transition |T^, and 
the other is that structure of phase spaces are understood in systems which are 
nearly integrable or fully chaotic. Models investigated in this paper consist of 
nonlinear oscillators with nearest neighbour interactions, and each oscillator is 
on a lattice point of the 3-dimensional simple cubic lattice. 

This paper is constructed as follows. We introduce five models in section 
1^. They are used to confirm the new universality of Lyapunov spectra appears 
in a wide class of Hamiltonian systems. We show Lyapunov spectra yielded 
by using random matrices with temporal 6- or exponential correlations in 
section ^. The Lyapunov spectra are straight even in systems being on the 
3-dimensional lattices when the systems are fully chaotic. In sections ^ and 
1^, we show that Lyapunov spectra for the five systems are not straight and 
have universality, and that the systems are neither fully chaotic nor nearly 
integrable. In section ^ we show that the new universality appears when 
quadratic terms of a potential function, U2, dominate higher terms, U4, namely 
U2/U4 takes large values. Section is devoted to summary and discussions. 
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2. Models 

We introduce five model Hamiltonians each of which represents a system being 
on the 3-dimensional simple cubic lattice with nearest neighbour interactions 
and periodic boundary condition. All Hamiltonians consist of kinetic and 
potential terms, 

H{q,p)=K{p) + U{q), (1) 
where the kinetic term is 

TV . 

K{P) = E :;pI (2) 

and is degrees of freedom, namely N = where L is the linear size of the 
lattice. 

One of the models is called XY model and expressed as follows 

Uxriq) = E [1 - cos(g. - g,)], qj G [0, 2tt), (3) 

<ij> 

where the summation J2<ij> takes over all the pairs of nearest neighbour lattice 
points i and j. 

The following three systems are expressed as 

1 ^ 

U{q)= E :^fe-?.)' + Eng.), (4) 
<ij> ^ j=i 

and they are distinguished by their potential functions V{q). In Double Well 
(DW) model 

VDwiq) = -\q' + \q\ (5) 



in Single Well (SW) model 



Vsw{q) = \q' + \q\ (6) 
and in Lorentzian (LO) model 

yL0{q) = (7) 
The last one has interactions of FPU-/3 type (3DFPU) 



UzDFPU — E 



k 1 

i^i.qi - qjf + ^{qi - qjY 



<i]> 

This model is used in section ^ to determine which term is dominant when the 
new universahty appears. 
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Numerical integrations of Hamiltonian equations of motion, 

_ dH{q,p) dpj_ _ dH{q,p) _ 
dt dp, ' dt dq, ' ^'^---'^^^ 

are performed with fourth order symplectic integrator with the fixed time slice 
At = 0.01. Accuracy of total energy is /^E / E ~ 0((At)*^) where and E 
are error and an initial value of total energy, respectively. 



3. Lyapunov spectra with random matrices 

As mentioned in the Introduction, Lyapunov spectra L{i/N) calculated with 
random matrices are straight in 1-dimensional chains. Lyapunov spectrum is 
a set of Lyapunov exponents {Aj} (i = 1,2, ■ ■ ■ , D), where D is dimension of 
phase space and the exponents are put in order as Aj > Aj+i. The summation of 
Aj up to n, Y17=i indicates linear instability of n-dimensional volume element 
Vn in phase space along a sample orbit. Namely, Vn diverges or converges as 

~exp[(Ai + A2 + --- + A„)t], (10) 

where t represents time. Each of Hamiltonian systems with degrees of 
freedom has 2N Lyapunov exponents which satisfy the following relations 
induced by symplectic properties 

A2iv-m = -A. it = l,2,---,N). (11) 

Hence we have only to observe the first half of Lyapunov spectrum when we 
consider Hamiltonian systems. Details of Lyapunov spectrum are reviewed in 



Ref. and references therein. 

The purpose of this section is to show that fully chaotic systems have 
straight Lyapunov spectra even in the 3-dimensional simple cubic lattice. We 
describe how we calculate Lyapunov spectra with random matrices, and then 
the Lyapunov spectra are shown. 

Lyapunov spectra indicate linear instability of a sample orbit, and they 
are calculated from linearized equations of motion, 

(12) 

where {Sq,6p) is a tangent vector, Ia? is the unit matrix of N x N, and (z,j) 
element of the matrix A{q) is 
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Here we used the form of our Hamiltonians, equations (|1|) and (0). Temporal 
evolution of the matrix A{q) is determined by temporal evolution of gj (t)'s, and 
we assume that qj{t)'s are independent random variables with S- or exponential 
correlations, namely 

Cij{t) oc 6ij6{t) or 6ije-'^\ (14) 

Cij{t) is the correlation function between qi{t) and qj{t), and it is defined as 

Cijit) = Aqi{t)Aq,{0) = lim ^ T dt'Aqi{t + t')Aq,{t'), (15) 

and 

Aqjit) = q,{t) - hm ^ f dtq,{t). (16) 

We show Lyapunov spectra in figure |^ which are calculated with random 
matrices having the 5-correlation, and we find that the Lyapunov spectra are 
straight as they are in 1-dimensional chains. We take random variables from 
the uniform distribution, and ranges of qj{t)^s are [— 7r,7r) for XY model, and 
[-3, 3] for DW and SW models. 



0.25 
0.2 

0.15 
0.1 

0.05 




(a) 



*V XY 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
i/N 



0.01 
0.008 
0.006 
0.004 
0.002 














Ox. 




SW ^v^v 


(b) 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
i/N 



(a) XY 



(b) DW and SW 



Figure 1. Lyapunov spectra yielded by using random matrices with S- 
correlation for XY, DW and SW models, (a) XY model, (b) DW «>) and 
SW (+) models. They are straight even in systems on 3-dimensional lattices 
as they are in 1-dimensional chains. Random variables follow uniform 
distributions. The straight solid lines arc guides for eyes. 



Next we change (5-correlation into exponential, and we set a = 0.4, 0.6, 0.8 
and 1.0 where a is the reciprocal number of correlation time of qj{t) (see 
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equation (]T^)). Figure ^ shows Lyapunov spectra obtained by using random 
matrices with exponential correlation, and the Lyapunov spectra are also 
straight in the region of 0.2 ^ i/N < 1. Consequently, we suppose finite 
correlation time does not affect the straightness of Lyapunov spectra. 
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i/N 

Figure 2. Lyapunov spectra yielded by using random matrices with 
exponential correlation for DW. a, the reciprocal number of correlation 
time, is 0.4,0.6,0.8 and 1.0 from lower to upper. Spectra show straight 
behaviour as figure || in the region of 0.2 ^ i/N < 1. The straight solid 
lines are guides for eyes. 




Random matrices with S- or exponential correlation yield straight 
Lyapunov spectra, and fully chaotic systems hence give the straightness even in 
the 3-dimensional simple cubic lattice. This result is used later to distinguish 
our systems evolved by Hamiltonian equations of motion from fully chaotic 
ones. 

4. New Universality of Lyapunov spectra 

In this section, we show that new universality of Lyapunov spectra exists in 
systems which are neither fully chaotic nor nearly integrable and which are in 
the thermodynamic limit (A^ —>■ oo) through giving the following four results, 
(i) The degrees of freedom = 4'^ is high enough to reach the thermodynamic 
limit for Lyapunov spectra, (ii) Forms of Lyapunov spectra are invariant with 
respect to energy in each of the four models which are XY, DW, SW and LO in 
a finite range of large i/N. (iii) The invariant forms of the four systems are in 
well agreement, and they are not approximated by the straight line. Here we 
conclude that new universality of Lyapunov spectra is found, (iv) Appearance 
of the new universality is not limited in nearly integrable systems. 
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Lyapunov spectra for XY model are shown in figure ^(a) in which points 
and dots represent that the degrees of freedom are N = 4^ and 10^ respectively. 
Numbers in the figure are values of energy density E/N. The Lyapunov spectra 
do not correspond with each other even though they have the same energy 
density. 
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Figure 3. Dependence on degrees of freedom N. Points and dots are 
Lyapunov spectra for A'^ = 4^^ and 10'^ respectively. Numbers in the figure 
represent values of energy density E/N . (a) Non-scaled, (b) Vertical axis 
is scaled for N = and the scale factor 7 = 1/1.1 for all values of energy. 
In each energy, scaled Lyapunov spectrum for = 4^ is in well agreement 
with one for N — 



Let US note that the form of Lyapunov spectrum L{i/N) is determined 
by ratios between Lyapunov exponents Aj rather than their absolute values. 
Accordingly, we may uniformly scale Lyapunov spectrum up or down from 
L{i/N) to '~fL{i/N), where 7 is arbitrarily picked for each spectrum. To 
multiple L{i/N) by 7 corresponds to change time scale from t into t/7. Scaled 
Lyapunov spectra are shown in figure §(b), and their forms holds regardless of 
degrees of freedom, where scale factors are 7 = 1/1.1 and 1.0 for each spectrum 
in = 4^ and 10'^ respectively. Since the thermodynamic limit of Lyapunov 
spectra [16||l7| exists, we suppose that the system reaches the thermodynamic 
limit even = 4^ concerning Lyapunov spectra. We therefore set A^ = 4^, 
and use scaled Lyapunov spectra without comment hereafter. 

Invariance of L{i/N) with respect to energy is shown in figure ^ for the four 
models. In each of the models, Lyapunov spectra are in well agreement among 
various values of energy in a range of large i/N m middle energy regime. Scale 
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factor 7 and the ranges oi i/N where the invariance appears are arranged in 
tables ^l| and |A2|, respectively. We remark that the invariance breaks or is 



strictly limited in a narrow range of i/N when energy is too high or low, that 
are E/N = 100 in XY and DW, and E/A^ = 1.0 in SW. 
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Figure 4. Lyapunov spectra for various values of energy in four 
Hamiltonian systems, (a) XY model, (b) DW model, (c) SW model, 
(d) LO model. Forms of Lyapunov spectra are the same in a range of large 
i/N in middle energy regime in each model. Scale factor 7 and the range 



of i/N giving the invariant form are arranged in tables Al and A2. 



Next we show the invariant forms of Lyapunov spectra of the four models 
together in figure ^. Scale factor 7 is arranged in table The four spectra 
in figure ^ coincide in a range of large i/N (0.4 ^ i/N < 1), while the straight 
solid line, which is obtained with random matrices, approximates them in a 
narrower range (0.7 ^ i/N < 1). That is, the invariant forms do not depend 
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on details of models and are not straight, and consequently, Lyapunov spectra 
of the four models have universality which is different from the one obtained 
in fully chaotic systems. 
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i/N 



Figure 5. Universal behaviour of scaled Lyapunov spectra of four models. 
The Lyapunov spectra coincide well among the four models in a range of 
large i/N, namely 0.4 < i/N < 1, and they are not approximated by the 
straight line, which approximates them only in the range of 0.7 H i/N < 1. 
The straight line is obtained with random matrices. Scale factor 7 is 
arranged in table A3. 



Moreover, we show that the new universality appears even when systems 
are not nearly integrable. We confirm that systems are not nearly integrable 
through showing that KAM tori, many of which survive in nearly integrable 
systems [jl8|, are not observed effectively. We use the fact that linear instability 



is suppressed around KAM tori since an orbit behaves like a regular one, 
while enhanced in chaotic sea. Accordingly, in nearly integrable systems, 
intermittency of local Lyapunov exponent \^°'^{n) occurs which is defined as 
follows. 

1 f{n+l)T 

A'r(n) = - / Ai(t)rft, (17) 

Ai(t) = ^iog|X(t)|, (18) 

where X{t) = {Sq, 6p) is a 2N dimensional tangent vector following linearized 
Hamiltonian equations of motion, equation (|12|). Figure ^ shows two time series 
of local Lyapunov exponent for E/N = 1.0 and 3.0 in XY model, both of which 
yield the new universal Lyapunov spectra. Here N = 10^. Intermittency is 
found for E/N = 1.0 but not for E/N = 3.0, and hence existence of KAM 
tori does not seem to be related to the new universality. Consequently, the 
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universality appears in systems with moderate strength of chaos which are 
between nearly integrable and fully chaotic. 
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Figure 6. Temporal evolutions of local Lyapunov exponent X^{"^{n) in 
XY model. N — 10"^. We set r — 10. Lower and upper curves take 
E/N — 1.0 and 3.0 respectively, both of which yield the new universal 
Lyapunov spectra. Intermittency is found for E/N ~ 1.0 but not for 
E/N = 3.0. 



5. 3DFPU model and quadratic interactions 



Table 1. Behaviour of Lyapunov spectra in 3DFPU model, k and E/N 
are coupling constant and energy density, respectively. The sign S means 
straight behaviour of the Lyapunov spectrum, while the sign C curved 
behaviour. The forms of curved Lyapunov spectra are in well agreement 
with the universal form. 
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To probe what conditions produce the universality, we study 3DFPU 
model with various values of energy density E/N and coupling constant k 
(see equation (^). Lyapunov spectra for 3DFPU model are shown in figure 
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Table 2. Ratios between time averages of the quadratic term U2 and the 
quartic terms U4 of potential function in 3DFPU model. We subtract 1.50 
from each of the ratios to make a threshold clear, namely values arranged 
in this tables are U2/U4 — 1.5. The positive values are found at the places 
where the sign C appears in table |l]. 
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1^ with a Lyapunov spectrum for XY model belonging in the universality, and 
they are classified into group C and group S. Group C includes the universal 
Lyapunov spectrum which is in well agreement with curved three spectra for 
3DFPU, and other spectra for 3DFPU are straight and belong to group S. 
The classification is arranged in table |T] with respect to k and E/N. Lyapunov 
spectra show curved forms when k is large and E/N is low, we hence conjecture 
that quadratic terms of potential function are dominant rather than quartic 
terms when the universality appears. 
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Figure 7. Lyapunov spectra in 3DFPU model. They are classified into 
group C and group S. Group C includes the universal spectrum which is 
in agreement with three spectra for 3DFPU, while group S consists of six 
spectra for 3DFPU which are straight. The straight line is a guide for eyes. 
The classification is arranged in table |l| with respect to k and E/N. Scale 
factor 7 shown in table A.4. 
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This conjecture is verified by taking ratios between time averages of U2 
and f/4, namely 1/2/114^, wliere 

U2=Y: % - q,)\ (19) 

<ii> ^ 

f/4= E (20) 

and 

Vn= dtVn{t), (n = 2,4). (21) 

Tlie ratios are arranged in table and they are large at the places where the 
sign C appears in table |l|. Consequently, the conjecture is verified. 

6. Summary and discussions 

In order to consider structure of phase spaces in Hamiltonian systems with 
moderate strength of chaos, we numerically investigated Lyapunov spectra 
{Aj} for five Hamiltonian systems with many degrees of freedom. We showed 
existence of universality of Lyapunov spectra, which is defined as Lyapunov 
spectra approximately take the same form regardless of energy and details of 
systems. The universality gives a curved form for Lyapunov spectra, and hence 
it is different from the one obtained in fully chaotic systems. 

A feature of the new universality is that it appears in a finite range of large 
ijN , where is degrees of freedom, and accordingly, properties depending on 
energy or details of models affect forms of Lyapunov spectra only in the range 
of small ijN where is out of the universality. In other words, we have only 
to focus on the range of small ijN when we are interested in such individual 
properties. 

We studied what conditions induce the universality to understand the cause 
of it. We showed that quadratic terms of a potential function dominate higher 
terms when the universality appears, harmonic motions are therefore regarded 
as the base of global motions in phase spaces. 

We geometrically interpret the harmonic motions and that the universality 
appears in a finite range of ijN as follows. Harmonic motions occur in high 
dimensional subspaces of phase spaces, which correspond to the finite range 
of ijN , and structure of phase space consists of chaotic sea and wrecks of 
n-dimensional tori, where n < N and n may change for each torus. Here n- 
dimensional torus means, roughly speaking, direct product of and {N — n)- 
dimensional instability, which is hyperbolic or complex. Note n = iV is KAM 
torus. 
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We need further analyses to confirm whether this interpretation is vahd 
or not, and to understand the origin of the new universality. We give two 
approaches which are geometrical and analytical. 

A geometrical approach uses resonance of instability along an orbit. 
Negative curvature of potential function induces positive Lyapunov exponent. 
On the other hand, even the curvature is always positive, resonance of 
instability along a sample orbit gets the largest Lyapunov exponent to be 
positive [T^. Here let us assume this resonance theory can be extended 



not only to the largest Lyapunov exponent but to all the exponents. If the 
resonances are yielded by harmonic motions in high dimensional subspace of 
phase space, then universality of Lyapunov spectra may be obtained because 
harmonic motions are independent of details of models. 

An analytical approach is to calculate decay rate spectrum of harmonic 
motions. In a dissipative system, behaviour of Lyapunov spectrum agrees with 
the decay rate spectrum of the linear fluctuation modes from the stationary 
solution [|10|. This suggests that to analyze decay rate spectrum is useful 



to understand behaviour of Lyapunov spectrum. Hamiltonian systems are 
regarded as dissipative systems when we observe only subspaces of phase 
spaces, and hence the decay rate spectrum must be useful because the new 
universality appears only in high dimensional subspaces of phase spaces. 

Our goal is to understand global structure of phase spaces in Hamiltonian 
systems with moderate strength of chaos. We must be allowed to come near the 
goal by understanding the cause of the new universahty, because ratios between 
Lyapunov exponents, which determine the form of the Lyapunov spectrum, 
seem to concern global structure of phase spaces. The new universality of 
Lyapunov spectra is an important clue to reveal the global structure. 
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Appendix A. Tables of scale factor 7 



Table Al. Scale factor 7 in figure ^. 
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2.0 


1.0 


1.4 


2.4 


3.0 


2.65 


3.0 


1.2 


3.0 


1/2.6 


1.8 


1.4 


5.0 


1.6 


5.0 


1.4 


4.0 


1/5.0 


2.2 


1.0 


7.0 


1.2 


7.0 


1.95 


5.0 


1/5.5 


100 


1/11 


9.0 


1.0 






100 


4.0 















Table A2. Rough estimations of ranges of i/N where universal behaviour 
of Lyapunov spectra appears in four models. 



Model XY 



DW 



SW 



LO 



i/N [0.4,1] [0.31,1] [0.16,1] [0.11,1] 



Table A3. Scale factor 7 in figure ^. 



Model XY DW SW LO 



E/N 

7 



6.0 2.0 3.0 3.0 
1.0 18 12.7 33 
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Table A4. Scale factor 7 in figure ||. 



k\E/N 

0.4 
0.7 
1.0 



1.0 2.0 3.0 



1.0 
1.7 
2.25 



0.75 
0.91 
1.4 



0.66 
0.74 
0.92 



